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We briefly discuss the state of the art on the anomalous dynamics of the Hamiltonian Mean Field 
model. We stress the important role of the initial conditions for understanding the microscopic 
nature of the intriguing metastable quasi stationary states observed in the model and the connections 
to Tsallis statistics and glassy dynamics. We also present new results on the existence of metastable 
states in the Kuramoto model and discuss the similarities with those found in the HMF model. The 
existence of metastability seem to be quite a common phenomenon in fully coupled systems, whose 
origin could be also interpreted as a dynamical mechanism preventing or hindering sinchronization. 
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I. INTRODUCTION 

The Hamiltonian Mean Field (HMF) model intensively 
studied in the last years can be considered a paradig- 
matic system for understanding the behavior of long- 
range interacting systems 0. In particular its strik- 
ing out-of-cquilibrium dynamics has raised much inter- 
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the anomalous dynamics connected with the existence of 
metastable quasi stationary states (QSS). Several claims 
have been advanced concerning a theoretical description 
of thisphenomenon within standard statistical mechan- 
ics 0, IE S • However these studies neglects the im- 
portant role of the initial conditions and the hierarchical 
fractal-like structures which are generated in the /x-space. 
In our opinion the situation is still much more complex 
and not completely clear although some progress has cer- 
tainly been made. In this very short paper we want 
summarize our point of view based on the most recent 
numerical simulations |14| which indicate a connection 
to Tsallis generalized statistics and glassy dynam- 
ics ^3|- We present also another interesting point of 
view, which could probably help in understanding this 
anomalous behaviour, exploiting an analogy with recent 
numerical results on metastability recently found in the 
Kuramoto model. In other words, metastability could 
likely be seen also as a kind of dynamical hindrance to 
synchronization. 



II. METASTABILITY IN THE HMF MODEL 

The HMF model describes a system of planar rotators 
or spin vectors ~s*j = (cos 6*;, sin 6i) with unitary mass and 
interacting through a long-range potential. The Hamil- 



tonian can be written as 
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being 9i and ^conjugate variables. The latter can also 
describe a system of fully interacting particles rotat- 
ing on the unit circle. Introducing the magnetization 

M = ji\Yli=i~^i\ as or der parameter, it is possible to 
solve the model exactly at equilibrium, where a second- 
order phase transition exists. The system passes from a 
low-energy condensed (ferromagnetic) phase with mag- 
netization M 7^ 0, to a high-energy one (paramagnetic), 
where spins are homogeneously oriented on the unit cir- 
cle and M = 0. The critical point is found at a tem- 
perature T c = 0.5, which corresponds to the critical en- 
ergy density U c — 0.75. The caloric curve is given by 
U = H/N= | + M 2 )Q. 

At variance with this equilibrium scenario, the out-of- 
equilibrium dynamics - explored by means of numeri- 
cal simulations - shows several anomalies in a special 
subcritical region of energy values and in particular for 
0.68 < U < U c . Starting from out-of-equilibrium ini- 
tial conditions, the system remains trapped in metasta- 
ble quasi-stationary states with vanishing magnetization 
at a temperature lower than the equilibrium one, until 
it slowly relaxes towards Boltzmann-Gibbs (BG) equi- 
librium. This transient QSS regime becomes stable if 
one takes the infinite size limit before the infinite time 
limit 0, Q ■ In this case the system stays forever in the 
M = state at the limitingtemperature Tqss- 
In several previous works @, E3] we have clearly shown 
that the microscopic nature of the anomalous QSS regime 
depends in a sensitive way on the choice of the initial 
conditions (IC). In this respect we focused on two main 
classes of IC (see Fig 1). 

(1) The first class of IC contains the so-called Ml IC, with 
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FIG. 1: Schematic figure illustrating the relaxation to BG equilibrium in the HMF model. According to the initial condition 
of the magnetization the nature of the QSS presents different features and the routes to thermalization are not the same. See 
text 



all the angles put equal to zero, i.e. M(0) = 1 and a wa- 
ter bag (uniform) distribution of momenta. In this class 
we can also consider those initial conditions with finite 
magnetization, M(0) > 0. When the system is started 
with these initial conditions an initial thermal explosion, 
with a following violent relaxation, occurs. After this 
stage, the system quickly freezes and remains trapped 
in long-living QSS with vanishing magnetization and a 
temperature smaller than the equilibrium one. In this 
case the system appears macroscopically quite homoge- 
neous, since averaging over all angles the magnetization 
tends to vanish with the size of the system. On the other 
hand from a microscopic point of view it is not so. In 
fact fractal-like structures are generated in the /i-space 
and many dynamical anomalies appears, among which 
anomalous diffusion, power-law decay of velocity corre- 
lations, aging, weak ergodicity breaking and dynamical 
frustration [2, 0- Such a behavior can be usefully de- 
scribed by the Tsallis' generalized formalism |l3.llq and 
interesting analogies with glassy systems do also exist 

(2) The second main class of out-of-equilibrium initial 
conditions indicated here as MO IC, see Fig.l, is charac- 
terized by a uniform distribution of both angles (M(0) — 



0) and momenta (water bag). This choice puts the system 
directly in the limiting long-living QSS with zero magne- 
tization and temperature Tqss- This is a particular case 
which corresponds to a stationary solution of the Vlasov 
equation [a, IE SB- The latter is stable in the thermody- 
namic limit but metastable for finite sizes of the system. 
After the QSS transient, the system relaxes very quickly 
(almost exponentially) to the BG equilibrium 3]. It is 
important to stress that although macroscopically also in 
this case we observe metastability as in the M1IC case, 
this time the initial state is spatially very homogeneous 
from the beginning and no violent relaxation occurs. The 
force acting on each particle (spin) is zero since t = and 
all the microscopic anomalies and correlations observed 
for initial conditions with finite magnetization are almost 
absent [jj. 

The very different microscopic behavior of the QSS 
regime corresponding to these two class of initial con- 
ditions strongly suggests a different nature of the macro- 
scopical metastability. In the MOIC case the Vlasov sta- 
bility argument 0,0 is rigorous due to perfect homogene- 
ity of the limiting temperature QSS. On the other hand 
in the M1IC case, microscopic fractal-like structures cre- 
ated by the sudden dynamical cooling of the system do 
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not correspond to the starting homogeneity hypothesis 
on which the Vlasov approach used in 0, ||| is based 
and create a much more complex situation [si El- In 
this latter case, the QSS metastability seems to be more 
related to the ergodicity-breaking phenomenon and to 
the trapping effect linked to the complex structure of the 
phase space region selected by the finite magnetization 
initial conditions. 

We summarize below the main numerical results on which 
our analysis is based. In refs.Q and elsewhere we studied 
the behavior of both velocity autocorrelation functions 
and anomalous diffusion in the QSS regime for U=0.69 
and for different initial conditions. A similar slow power- 
law decay of velocity correlations for 0.4 < M(0) < 1 
was found, while for M(0) = 0.2, and even more for 
M(0) = 0, the decay is much faster. By fitting these 
relaxation curves by means of the Tsallis' q-exponential 
function it was possible to characterize in a quantita- 
tive way the dynamics originatined by these different ini- 
tial conditions. In fact we found a value q = 1.5 for 
M(0) > 0.4, while we got q = 1.2 and q = 1.1 for 
M(0) = 0.2 and for M(0) = respectively - notice that 
q = 1 corresponds to the usual exponential decay. On the 
other hand, by studying the mean square displacement 
of angles (which typically scales as a 2 (t) ~ f), again for 
U=0.69 and as a function of the initial conditions, it was 
found superdiffusive behavior for 0.4 < M(0) < 1 with 
an exponent 7 = 1.4 — 1.5. At variance, for M(0)=0 we 
obtained 7 = 1.2, a value that tends to 7 ~ 1 (normal 
diffusion) increasing the size of the system. 
In refs. 0] it was found a connection between the velocity 
correlations decay and the diffusive behavior of metasta- 
ble states based on the simple relationship 7 = 2/(3 — q), 
which links the diffusion exponent 7 with the entropic 
index q. This formula was derived theoretically by Tsal- 
lis and Bukman within a generalized Fokker-Planck ap- 
proach based on q-statistics . For various initial con- 
ditions, ranging from M(0)=1 to M(0)=0 and different 
sizes at U = 0.69, it was found that this formula is cor- 
rect. This latter holds also for a generalized version of 
the model, the a — XY model, where the long-range char- 
acter of the interaction is modulated through a term r~ a 
|12| . In such a way one can predict the value of q from 
the knowledge of 7 and vice- versa. 

Another framework that seems to be very promising in 
order to shed light on the microscopic nature of the QSS 
is that one of glassy dynamics. The first link between 
QSS metastability and glassy behavior was suggested by 
the discovery of the aging phenomenon in autocorrela- 
tion functions 0, 0] and by the dynamical frustration 
scenario observed in the QSS regime obtained for M1IC. 
In this case clusters appear and disappear on the unit 
circle, competing one with each other in trapping more 
and more particles and generating a dynamically frus- 
trated situation typical of glassy systems. In refs.[l3j 
it was introduced a new order parameter for the HMF 
model, the polarization, in order to measure the degree 
of freezing of the rotators. By means of this new pa- 



rameter one can discriminate between M1IC and M0IC 
metastability, showing that only starting from the first 
class of initial conditions the QSS regime could be con- 
sidered (in the thermodynamic limit) as a sort of glassy- 
phase for the HMF model. More recently we introduced 
an effective spin-glass Hamiltonian in order to study this 
anomalous glassy-like dynamics from an analitical point 
of view. By means of the Replica formalism, it was found 
a self-consistent equation for the glassy order parameter. 
The latter is able to predict, in a restricted energy re- 
gion below the phase transition, the polarization values 
obtained with molecular dynamics simulations |l3j . 
All these results indicate a deep dynamical nature of QSS 
metastability. In the next section we present new results 
concerning metastability in the Kuramoto model which 
shows a certain analogy with those of the HMF model. 

III. METASTABILITY IN THE KURAMOTO 
MODEL 

One of the the simplest and most successful models for 
synchronization is the famous Kuramoto model of cou- 
pled oscillators 0, . It is simple enough to be an- 
alytically solvable, still retaining the basic principles to 
produce a rich variety of dynamical regimes and synchro- 
nization patterns. The dynamics of the model is given 
by 

K N 

0/(i)= Wi + -^sm(0 i -0 i ) i = l,...,N (2) 

where 9i(t) is the phase (angle) of the «th oscillator at 
time t (— 7r < 9i(t) < 7r), while u>i is its intrinsic frequency 
randomly drawn from a symmetric, unimodal distribu- 
tion g(uj) with a first moment loq (typically a Gaussian 
distribution or a uniform one) . These natural frequencies 
u>i are time- independent. The sum in the above equa- 
tion runs over all the oscillators so that the model is an 
example of a globally coupled system. The parameter 
K > measures the coupling strength: for small values 
of K , each oscillator tends to run independently with its 
own frequency, while for large values of K, the coupling 
tends to synchronize (in phase and frequency) the os- 
cillator with all the others. Kuramoto showed that the 
model, despite the difference in the natural frequencies 
of the oscillators, exhibits a spontaneous transition from 
incoherence to collective synchronization, as the coupling 
strength is increased beyond a certain critical threshold 

In Fig. 2 we plot the asymptotic order parameter r = 
e l6j versus the coupling K. In the simulations 
we used a uniform distribution for the natural frequen- 
cies and a homogeneous initial conditions for the phases, 
i.e. r(0) = 0. The figure shows a sharp transition from 
a homogenous non-synchronized state to a synchronized 
regime above K c ^2.6. We considered N — 10000 and 
average over 10 runs. Also in this case one can observe 
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FIG. 2: For the Kuramoto model we plot the asymptotic order parameter r as a function of the coupling strenght K. We 
considered N=10000 and an average over 10 runs. The grey area indicates the region where we have found the presence of 
metastable states. See text and next figure. 




FIG. 3: For the Kuramoto model, we show the metastable states found for values of K taken inside the grey area of the previous 
figure. The number of oscillator is N=10000 and an average over 10 runs was considered. 



the emergence of metastable states just below the criti- 
cal point. The grey area indicates the region where this 
phenomenon has been found. In Fig. 3 we show the 
time evolution of the order parameter for different values 
of the coupling K in the metastability region. Also in 
this case the lifetime of the QSS seems to depend on the 
size of the system in close analogy with the metastable 
states found in the HMF model. A detailed study of this 
behavior is in progress and will be reported elsewhere. 

Despite the non hamiltonian character of eq.(3), there 
is a formal link between the Kuramoto and the HMF 



model. Actually, both models can be considered as par- 
ticular cases of the following general equation describing 
a generic pendulum driven by a constant torque T: 

A6" + B9' + Csin{6) = T (3) 

where B is a viscous damping constant while A and C are 
constants depending on mass, lenght and force. In fact in 
the conservative case {B = T = 0), from eq.© we recover 
the equations of motion of HMF model rewritten using 
the order parameter M, with the global angle $ = and 
considering each spin (pendulum) as moving in a mean- 
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field potential determined by the istantaneous positions 
of all the other spins (pendula) : 



0" = -Msin(0i) 



i = 1, 



,N. 



(4) 



On the other hand, if we consider the overdamped limit 
(when A is negligible with respect to the friction coeffi- 
cient B), Eq.® reduces to the Kuramoto equation, writ- 
ten also in this case using the order parameter r with the 
the global phase equal to zero: 



6i = u>i — K r sin(0i) 



i=l,...,N. 



(5) 



In the latter case the natural frequencies u>i play the role 
of the constant torque T. 

These formal similarities provide further support to the 
possibility of a common explanation of metastability in 
terms of a dynamical hindrance of synchronization. 



of metastability for the HMF model and the fundamen- 
tal role played by the initial conditions, which are able 
in some cases to produce a non homogeneous dynam- 
ics hardly tractable within a Vlasov approach. This has 
lead us to follow several approaches which differ from 
those used in refs. [U IS 0> 0] • Tsallis statistics, glassy 
dynamics and hindrance of synchronization provide very 
interesting perspectives and promising routes to be fur- 
ther explored in the future for a better understanding of 
the dynamical origin of metastability and complex be- 
haviour in fully coupled systems. 



Acknowledgments 



IV. CONCLUSIONS 

Concluding this short review, we want to stress again 
the importance of the dynamics to understand the origin 
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